Draft version February 26, 2013 

Preprint typeset using I^T^X style cmulatcapj v. 5/2/11 



(N 

o 

(N 



o 
u 

6 



> 

ON 

o 
en 

en 

o 

CN 



13 



ULTRALUMINOUS STAR-FORMING GALAXIES AND EXTREMELY LUMINOUS WARM MOLECULAR 

HYDROGEN EMISSION AT Z = 2.16 
IN THE PKS 1138-26 RADIO GALAXY PROTOCLUSTER 

P. Ogle 1 , J. E. Davies 1 , P. N. Appleton 2 , B. Bertincourt 3 , N. Seymour 4 , G. Helou 1 

1 IPAC, California Institute of Technology, Mail Code 220-6, Pasadena, CA 91125, ogle@ipac.caltech.edu 
2 NHSC, California Institute of Technology, Mail Code 220-6, Pasadena, CA 91125 

Columbia University Astronomy Laboratory and 
4 CSIRO Astronomy & Space Science, P.O. Box 76, Epping, NSW 1710, Australia 
Draft version February 26, 2013 

ABSTRACT 

A deep Spitzer Infrared Spectrograph (IRS) map of the PKS 1138-26 galaxy protocluster reveals 
ultraluminous PAH emission from obscured star formation in three protocluster galaxies, including 
Ha-emitter (HAE) 229, HAE 131, and the central Spiderweb Galaxy. Star formation rates of ~ 
500 — 1100A/ Q yr _1 are estimated from the 7.7/im PAH feature. At such prodigious formation rates, 
the galaxy stellar masses will double in 0.6-1.1 Gyr. We are viewing the peak epoch of star formation 
for these protocluster galaxies. However, it appears that extinction of Ha is much greater (up to 
a factor of 40) in the two ULIRG HAEs compared to the Spiderweb. This may be attributed to 
different spatial distributions of star formation-nuclear star formation in the HAEs versus extended 
star formation in accreting satellite galaxies in the Spiderweb. We find extremely luminous mid- 
IR rotational line emission from warm molecular hydrogen in the Spiderweb Galaxy, with L(H2 0-0 
S(3))= 1.4 x 10 44 erg s _1 (3.7 x 10 10 i©), ~ 20 times more luminous than any previously known 
H2 emission galaxy (MOHEG). Depending on temperature, this corresponds to a very large mass of 
> 9 x 10 6 - 2 x 10 9 A/ Q of T > 300 K molecular gas, which may be heated by the PKS 1138-26 radio 
jet, acting to quench nuclear star formation. There is > 8 times more warm H2 at these temperatures 
in the Spiderweb than what has been seen in low-redshift {z < 0.2) radio galaxies, indicating that the 
Spiderweb may have a larger reservoir of molecular gas than more evolved radio galaxies. This is the 
highest redshift galaxy yet in which warm molecular hydrogen has been directly detected. 
Subject headings: Galaxies:high-redshift, Galaxiesdormation, Galaxies:clusters:individual:PKS 1138- 
26 , Galaxies : individual: Spiderweb 



1. INTRODUCTION 

Galaxy clusters are the largest virialized structures 
in the universe. Massive, X-ray se lected clusters are 
well established at redshift z = 1 (Finogu cnov et al. I 
l2007f ). and their brightest galaxies hav e colors consistent 
with passively evolving ellipticals (e.g.. lMcCarthv et al. I 
l2007h . The progenitors of such clusters must have formed 
at an earlier epoch, around the highest density peaks 
in the cosmos. Such protoclu sters have been found 
around z = 2 — 5 radio galaxies (lPentericci et al." 200 
Venemans et al. | [2007; Hat ch et al. 11201 lHKuiper et al 



201lUTanaka et al. 11201 II) . which appear to be construe 



tion zones of cluster ellipticals and their su permassive 
black holes (|Lacv. Bunker, fc Ridgway1l2000t l. 

The baryonic mass of z < 1 galaxy clusters is dom- 
inated by primordial gas that is virialized and heated 
to 10 7 K in the potential wells of d ark matter halos 
(jAndreon I |20lct iGiodini et al. 1 I2009D . A reduction of 
the gas fraction relative to the cosmic baryon fraction can 
be explained by a combination of st ar formation, super- 
nova feedback, and AGN feedback (jBode et al. 1 120091 ) . 
There is evidence that z — 2 protoclusters are X-ray de- 
ficient ( Lx < 10 44 erg s" 1 ) com pared to clusters at lower 
redshift (|Overzier et al. |[2005T) , supporting the idea that 
they are in the process of formation and not yet virial- 
ized. 

Chandra X-ray imaging of the Spiderweb radio galaxy 



(PKS 1138-26, z — 2.156) shows e xtended X-ray emis - 
sion aligned with the radio lobes (|Carilli et al. I [2002) . 
which indicates that the X-ray emission may come from 
gas that is shock-heated by the radio jet rather than viri- 
alized in the cluster potential. Large polarization ro- 
tation measures are found for PKS 1138-26 and other 
high-redshift radio galaxy lobes, indicating the pres- 
ence of large column densities of ionized, ma gnetized 
gas (|Carilli et al. I Il997t lAthreva et al~l 119981 ). This 

,as may have it s origin in radio jet-driven outflows 

Nesvadba et"aT~ll205g ). 

Large quantities of shock-heated warm molecular gas 
are foun d by Spitzer in 30% o f 3C radio galaxies a t 
z < 0.2 (|Qgle et al. I [20101 [20071 IGuillard et al. 1 120121) . 
The heretofore most luminous and distant H2 emis- 
sion is seen in the z = 0.29 Zw 3146 brightest clus- 
ter galaxy, with a single H2 line luminosity L(H2 0-0 
S(3))= 6.1 x 10 42 erg s" 1 (Egami et al. I [20061) . In ra- 
dio galaxies, the H2 is most likely powered by kinetic 
energy from the relativistic jet that is dissipated in the 
host galaxy interstellar medium (ISM). The large ratio 
of H2 to PAH luminosity observed in these and other 
H2 emission galaxies (MOHEGs) can not be generated 
in stellar photodissociation regions (PDRs) . X-ray heat- 
ing by relatively low-luminosity AGNs is also insufficient 
to power the obs erved H2 rotational line luminosities 
(jQgle et al. Il2010l) . Star forma tion appears to be su 



pressed in radio-loud MOHEGs (jNesvadba et al 
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falling a factor of ~ 50 below the Kennicutt-Schmidt rela- 
tion bet ween specific star -formation rate and gas surface 
density (jKennicutt II1998D . 

A spectacu lar protocluster su rrounds the Spiderweb 
radio galaxy ([Kurk et al. I [2000). As many as 121 can- 
didate cluster members within a 4.5' (2.2 Mpc) ra- 
dius have been identified as possible Lya emitters (37 
LAE), Ha emitters (40 HAE), or extre mely red ob- 
jects (44 E ROs) via filter photometry (|Kurk et al. I 
2000, 2004a). A large overdensity of sources with red 
IRAC [3. 6] -[4. 5] colors or bright 24 / nm emission are 
found within a radius of ~ 1 — 2' (Galametz et al. I 
120121: iMavoet al.~ll2012D . At least 24 cluster members 
have been spectroscopically detected in Lya (15) or Ha 
emission (9 ) and have reds hifts close to that of the 
Spiderweb (IPentericci et al. I [20 00; Kur k et al. I I2004ri 
IVenemans et al. H2007t IDohertv et al. 112010ft . The clus- 
ter consists of two merging sub-clumps, with velocity 
dispersions of 520 and 280 km s _1 for the associated 
LAEs and a com bined virial mass of ~ 1 x 10 14 M Q 
(|Pentericci et al. 1 12000). Several HAEs are bound to the 
Spiderweb Galaxy with a velocity dispersio n of 360 km 
s -1 , in dicating a virial mass of 8 x lO 13 M0( |Kurk et al. I 
I2004bft . 

Rest UV to optical continuum imaging of the clus- 
ter core with HST shows multiple satellite galaxies with 
peculiar mor phologies within 150 kpc of the Spider- 
web Galaxy (IPentericci et al. I [19981 IMilev et ai~ll2006t 
IHatch et al. I I2009D . An enormous Lya halo imaged 
by the VLT surrounds the galaxy out to a radius of 
> 40 kpc, indicating a large reserv oir of atomic gas 
(IKurk et al. l[200llHatch et al. 112008ft . 

We made a deep Spitzer infrared spectroscopic map 
of the central region of the PKS 1138-26 protocluster, 
centered on the Spiderweb radio galaxy, in order to ex- 
plore mid-IR indicators of star formation and shocked 
gas. Here we present mid-IR spectra of the Spiderweb 
and two other protocluster galaxies which are bright in 
the MIPS 24 (im band. Spectra of the remaining sources 
will be presented in a future paper. Distance-dependent 
quantities in this paper are calculated using cosmological 
parameters H = 70 km s -1 , il\ = 0.7, and Q m = 0.3. 

2. OBSERVATIONS 

We used the Spitzer Infrared Spectrograph 
(|Houck et al. I I2004H to map a 3.2' x 1.1' strip of 
the PKS 1138-26 protocluster (Fig. 1). The LL1 (20.5- 
38.5 /zm) low-resolution slit (168" x 10'.'7) was stepped 
perpendicular to the slit in 11 x 5'.'0 steps and parallel 
to the slit in 5 x 5'.'0 arcsecond steps to create a map 
centered on the Spiderweb radio galaxy. Exposures were 
2x120 s per pointing. The map was repeated 8 times 
over 5 days, giving a total exposure time of 29.3 hr and a 
mean exposure time per pixel of 4.2 hrs. Rotation of the 
spacecraft and spectrograph slit (PA= 124.6 — 126.3°) 
over this time period improved the spatial sampling of 
the map. Parallel data were also recorded in the LL2 
(14-21 (im) slit, which will be reported elsewhere. 

Spectra (Fig. 2) were reduced from SSC Pipeline 18.18 
basic calibrated data (BCDs) using pipeline uncertainties 
and custom background subtraction. The background 
was estimated by stacking the BCD frames in a time se- 
ries, masking sources, and fitting each pixel time series 
with a robust, low-order polynomial. The LL1 spectral 



cube was built using CUBISM ([Smith et al. Il2007aft . us- 
ing a spatial oversampling factor of 3, and cleaned us- 
ing automatic 10c and interactive bad-pixel identifica- 
tion and removal. A noise cube was constructed from 
the difference of the first and second exposure for each 
pointing, and its pixel-to-pixel standard deviation was 
found to agree closely with the pipeline-derived uncer- 
tainties as a function of wavelength. For most purposes, 
we confidently use the pipeline uncertainties in our data 
analysis. However, these uncertainties appear to signifi- 
cantly underestimate the noise at observed wavelengths 
> 33 /xm (rest frame 10.4 /an). We attribute this to sys- 
tematically deviant pixels that we were unable to iden- 
tify and remove in this noisier wavelength range. Feature 
fluxes measured in this wavelength range, such as those 
of [S iv] and the 11.3 /im PAH, should be interpreted 
with caution. 

Sources were identified and their positions measured 
in the collapsed IRS cube, weighted by the MIPS 24 
[im band spectral response (Fig. 1). Spectral extrac- 
tions were performed using a 3D optimal extraction al- 
gorithm that we adapted fr om SPICE 2D optimal ex- 
traction ([Narron et al. 1 2007; Home Ill986[ ). We used a 
wavelength-smoothed and interpolated 3D PSF template 
derived from the Spiderweb itself for optimal extraction. 
The flux from this source is dominated by the spatially 
unresolved AGN. Spectra were optimally extracted (us- 
ing S/N weighting) within an aperture that expands lin- 
early with wavelength and has a diameter of 8 pixels (2.7 
native LL1 pixels, or 14") at 27/im. Such a small aper- 
ture is utilized to minimize contamination from nearby 
sources. A wavelength-independent flux calibration fac- 
tor equal to the pixel solid angle (5.624 x 10~ n sr/pixel) 
was derived from the ratio of SPICE over CUBISM ex- 
tractions of the Spiderweb, in BCDs where it was cen- 
tered on the slit. We applied a correction for detector 
fringeing to the Spiderweb spectrum, reducing ~ 3% os- 
cillations in the observed 21-25 /im region down to the 
1% level. 

We used PAHFIT ([Smith et al. ll2007bD to obtain a 
least-squares fit to the continuum, emission lines, and 
PAH features. The continuum was fit with the sum of 
blackbody spectra covering a range of dust temperatures. 
No extinction was included in the fits since it did not im- 
prove them. The emission lines were fit with unresolved 
Gaussian line profiles. In the case of emission line blends, 
the fluxes were summed to give the total flux of the blend, 
rather than report individual line fluxes. The PAH fea- 
tures are fit by Drude profiles, with rest wavelengths and 
widths fixed to values appropriate for star-forming low- 
redshift galaxies, but with relative intensities allowed to 
vary. While it is not clear a priori that the mean PAH 
composition and feature shapes should be the same at 
high redshift, this provides a good fit for most features, 
except for the 8.6 /im PAH which appears to be shifted 
to shorter wavelength (§3). 

3. RESULTS 

We present the stacked IRS cube in Figure 1, weighted 
by the MIPS 24 /im band spectral response. We detect 13 
sources with peak flux > 0.2 mJy that yield high-quality 
spectra with peak S/N> 10 per 0.18 /im bin in the ob- 
served 21-35 /im band. In this paper, we concentrate on 
the 3 brightest 24 /im sources matched to spectroscopi- 
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cally confirmed cluster members-Spiderweb, PKS 1138- 
26:[KPR2004] HAE 0229, and PKS 1138-26: [KPR2004] 
HAE 0131 (hereafter HAE 229 and HAE 131). 

The IRS spectrum of the Spiderweb is dominated by 
hot dust emission from the central AGN (Fig. 2a, b). 
A silicate emission bump observed at 9.3-11/im (rest) 
is fit by a 120 K blackbody multiplied by the silicate 
emissivity function. We detect PAH emission at 7.7 /im, 

8.6 /jm, and 11.2 /jm. The 11.33 /im PAH feature falls 
in a noisy region of the spectrum, so it is formally not 
detected. However, the majority of the reported 11.23 
/im flux (Table 1) is likely contributed by the blue wing 
of the 11.33 /im PAH feature. 

We detect a number unresolved emission lines (Table 2) 
in the spectrum of the Spiderweb, including [Ne Vl], [Ar 
in], [S IV], and H 2 0-0 S(3). The H 2 0-0 S(4) line may be 
detected, but could be confused with substructure near 
the peak of the 7.7 /im PAH feature. A blend of [Ar n] 
and H 2 0-0 S(5) is detected at 6.9 /im. We confirm the 
H2 0-0 S(3) line detection with a continuum-subtracted 
map of the line at z = 2.156 (Fig. 1, bottom). This 
line is only detected at the location of the Spiderweb 
Galaxy, with a peak significance of > 6<r in the map and 
an integrated line detection of 14a in the spectrum. 

HAE 229 and HAE 131 both have PAH-rich spectra 
(Figs. 2c,d, and Table 1). We detect PAH blends at 

7.7 /im, 8.6 /im, and 11.3 /im (rest) in both galaxies. 
The [Ar 11] line is detected in HAE 229 and marginally 
detected in HAE 131. The [Ne vi] and H 2 0-0 S(4) lines 
may be detected in HAE 131, however it is difficult to 
clearly distinguish these lines from possible substructure 
in the 7.7 /im PAH feature. While the H 2 0-0 S(3) line 
is formally detected in HAE 131 (0.11 ±0.03 x 10" 14 erg 
s _1 cm -2 ), the peak of the line is not at the expected 
wavelength and may be contaminated by a noise spike, 
so we treat this as an upper limit. This line is clearly not 
detected in the continuum-subtracted H 2 S(3) map (Fig. 
1). Hot dust continuum underlies the PAH spectrum in 
HAE 229 and weaker continuum emission is also visible 
at the short wavelength end of the HAE 131 spectrum. 

Shifting the wavelength of the 8.6 /im PAH feature 
from the standard wavelength of 8.61 /im to 8.57 /im 
(corresponding to a Doppler shift of -1400 km s" 1 ) gave 
a better fit to the spectra of all three sources. We be- 
lieve this indicates a shift in the intrinsic wavelength 
rather than a Doppler shift because the emission lines 
and 7.7 and 11.3 /im PAH features have wavelengths con- 
sistent with the optical-NIR emission line redshifts. The 
8.6/7.7 micron PAH ratios (0.21-0.38) seen in the three 
galaxies are somewhat greater than the range of PAH ra- 
tios (0.11-0.21) seen i n low-redshift star-forming galaxies 
([Smith, et al. Il2007bl ). These differences in 8.6 /im PAH 
wavelength and PAH ratios may indicate a difference in 
chemical composition or excitation of the PAH molecules 
compared to low-redshift galaxies. 

4. DISCUSSION 

4.1. Spiderweb Galaxy 

The Spiderweb Galaxy is surrounded by a large num- 
ber of smaller satellite galaxies, seen in HST ACS images 
(iPentericci et al. 1119981; iMilev et ai~ll2006t IHatch et al. I 
l2009f) . Eleven cluster members or candidates have been 
identified within our 14" diameter spectral extraction 



aperture, according to NED (Fig. 1). These include the 
Spiderweb, X-ray source CXOGBA Jl 14047.9-262906, 4 
LAE candidates, 3 HAE candidates, and 2 EROs. Ex- 
cept for the Spiderweb, none have spectroscopic redshifts 
available. From the strong hot dust emission and low 
equivalent widths of the PAH features in the Spiderweb 
spectrum, it is clear that most of the observed 24 /im 
emission comes from the AGN. The strong [Ne vi] and 
[S iv] emission lines are additional indicators of a domi- 
nant AGN contribution to the spectrum. 

The Spiderweb has a 7.7/im PAH luminosity of 1.7 x 
10 45 erg s" 1 (4.4 x lO 11 ^, Table 3). I n com parison, the 
brightest ULIRG in the lElbaz et all (|20lH) sample has 
a PAH-dominated 8 /im luminosity of 2 x lO 11 ^©, star- 
forming IR (SFIR) luminosity of L S fir ~ 3.3 x 10 12 L Q , 
and a star formation rate of SFR = 500M Q yr" 1 . The 
corresponding conversions among 8 /im luminosity, total 
star-forming IR luminosity, and star formation rate are 
L S fir/L$ = 17 and i 8 /SFR= 4 x IO^qM" 1 yr. Lack- 
ing rest-frame far-IR fluxes, we adopt these as represen- 
tative conversion factors for the ULIRGs in our sample. 
This gives an approximate PAH-based star-forming FIR 
luminosity for the Spiderweb of Lsfir ~ 7.4 x 10 12 L© 
and a star-formation rate of ~ IIOOMq yr -1 . Equally 
large PAH luminosities and PAH-based star-formation 
rates have be en reported in the z = 3.53 radio galaxy 6C 
J1908+7220 (ISevmour et all 12003, 4OOOM yr" 1 ) the 
Cloverleaf Quasar (|Lutz et al. 1120071 ~ IOO OMq yr" 1 ) 
and other high-redshift mm-bright QSOs ([Lutz et al. I 
l200l . 

We estimate a similar star formation rate of ~ 1200Af Q 
yr -1 from the extended H a luminosity of the Spider- 
web (|Nesvadba et al. 1120061). Here we use SFR(Ha , M q 
yr" 1 )= 3.1 x 10" 8 i(Ha) ([Kennicutt et al. 1119941 1 and 
assume no extinction. The extended Ha luminosity was 
measur ed after subtracting th e point-like AGN contri- 
bution ([Nesvadba et al. II2006T) . but the AGN might still 
contribute to the extended Ha emission via photoioniza- 
tion or radio-jet driven shocks. 

A much lower star formation rate of 57 ± 8 M Q yr" 1 
is estimated f rom the diffuse 32-kpc radius UV halo 
(|Hatchet al. 1 120081 ). This rate increases to 142 M Q 
yr" 1 correcting for screen extinction (E(B-V)~ 0.1 from 
rest UV color), and 325 M Q yr" 1 if all of the diffuse 
and compact UV s ources are included in the estimate 
([Hatchet al. 1 120081 ) . This is still far below the PAH- 
based star-formation rate, indicating that the assump- 
tion of a uniform dust screen is too simplistic. Instead, 
much of the star formation in the Spiderweb must be 
highly obscured (> 2.2 mag) at rest-UV wavelengths. 

The Spiderweb is the brightest and most massive 
galaxy in the PKS 1138-26 cluster, with a stellar mass 
of 1.1 x 10 12 Mf7) estimated from the SEP of the centra l 
radio galaxy ([Hatch et al. 1120091 : ISevmour et al. 112007ft . 
Adding the masses of the satellite galaxies estimated by 
IHatch et al. I (|2009f ) increases the total stellar mass to 
1.2 x 10 12 M Q within our 7" (60 kpc) radius spectral ex- 
traction aperture. The stellar mass doubling time scale 
at the PAH-inferred star formation rate is therefore ~ 1.1 
Gyr. It appears that the Spiderweb is in the process 
of building the bulk of its stellar mass at the observed 
epoch. The Spiderweb may gain additional stellar mass 
(a factor of 1.1-2.2) by accreting known satellite galaxies 
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within a 150 kpc radius (|Hatch et al. II2009T) . 

We measure an H 2 0-0 S(3) luminosity of 1.4 x 10 44 
erg s" 1 (3.7 x 10 10 i Q ). We note that the H 2 0-0 S(2), 
S(l) and S(0) lines may be similarly luminous, but are 
redshifted beyond the longest wavelength detectable with 
the IRS. Taken at face value, the H 2 0-0 S(3) and S(5) 
line fluxes fit by PAHFIT would indicate a temperature 
of ~ 920 K and an H 2 mass of ~ 1 x 10 7 M©, assuming 
a single-temperature gas in equilibrium. However, the 
H 2 0-0 S(5) line is blended with [Ar n]. We estimate 
a more conservative upper limit on the H 2 temperature 
and lower limit on the H 2 mass by assuming that the [Ar 
ii] + H 2 0-0 S(5) emission line fit gives an upper limit on 
the H 2 0-0 S(5) line flux. This yields an H 2 temperature 
(for gas contributing to the S(3) and S(5) lines) of 300 
K< T < 1200 K and a warm molecular hydrogen mass 
of > 9 x 10 6 - 1 x 10 9 M @ (Fig. 3). 

In the same temperature range, 3C radio galaxies at 
z < 0.22 ((Ogle et al. 1120 lOh have warm H 2 masses more 
than 8 times smaller than the Spiderweb (Fig. 3). In 
the lowest temperature range probed by Spitzer, their 
H 2 0-0 S(0)-S(2) lines indicate warm (T = 100 - 300 K) 
H 2 masses > 8 x 10 6 - 2 x 10 10 M Q . It is likely that 
the Spiderweb has a much larger total warm H 2 mass 
at T = 100 — 300 K than what we have detected with 
the 0-0 S(3) and S(5) lines at T > 300 K. If the ratio of 
warm to total H 2 mass is similar to that of low-redshift 
radio galaxies, then the extremely luminous H 2 emission 
from the Spiderweb could be explained by an extremely 
large molecular gas mass of ~ 2 x 10 11 M Q . This would 
give a reasonable H 2 to stellar mass ratio of ~ 0.2. A 
similarly large molecular gas mass of 1.4 x 10 11 Mq has 
been inferred from CO (1-0) emission fr om the z = 2.8 
lense d quasar SMM J041327.2+102743 ()Riechers et al. I 
l20ll . 

We measure a lower limit (from S(3) only) for the ratio 
of the summed first four H 2 0-0 pure rotational lines to 
7.7 a™ PAH of i(H 2 0-0 S(0)-S(3))/L(PAH7.7)> 0.08, 
more than a factor of 2 greater than the empirical and 
theoretical ratio of < 0.04 for UV heating in a PDR. This 
makes the Spid erweb Galaxy the most l uminous known 
radio MOHEG (IQgle et al. 1120071 l20loT: IGuillard et aT~l 
120121 ). a class of galaxies where the molecular gas is 
heated by dissipation of mechanical energy. By anal- 
ogy to low-redshift radio MOHEGs, the presence of a 
large mass of warm molecular hydrogen indicates that 
a significant fraction of the molecular gas content of 
the Spiderweb has been shock-heate d by the radio jet , 
in turn suppressing st ar formation (IQgle et al. 120101 : 
iNesvadba et al. 1 12010( 1 . iNesvadba et al. I ((2006( 1 esti- 
mate a radio jet mechanical power, based on the radio lu- 
minosity of PKS 1138-26, of - 10 46 erg s" 1 . Only - 1% 
of this would be enough to power the observed H 2 0-0 
S(3) line. On the other hand, the absorption-corrected 
(intrinsic A H = 2.6 x 10 22 cm" 2 ) 2-10 ke V X-ray lumi- 
nosit y of the AGN is only 4 x 10 45 erg s^" 1 (|CariHi et al. I 
2002), which is insufficient to power the H 2 emission. 
The observed ratio L(H 2 0-0 S(0)-S(3))/L x (2-10 keV)> 
0.035 is more than 5 times greater than the theoreti- 
cal maximum ratio for AGN X-ray hea ting of < 0.007 
((Ogle et al. I I2010t IGuillard et al. I [2011 . 

From our Spitzer spectrum, the monochromatic mid- 
IR luminosity of the dust-obscured Spiderweb AGN is 



uL v (U (jaa)= 1.9 x 10 46 erg s _1 . Using an IR/UV flux 
conversion factor of /uv = 0.36 and bolometric IR cor- 
rection factor of /boi = 1.7 for the AGN dusty torus, 
both derived from z < 1 radio galaxies and quasars 
((Ogle et al. 1 12006( 1 . we estimate an AGN bolometric lu- 
minosity of L b oi = tfho\/ fw)vL v (V\. ^m)= 9.5 x 10 46 
erg s _1 , ~ 3 times the total IR luminosity from star 
formation. This is a factor of ~ 2 lar ger than th e bolo - 
metric luminosity estimate of INesvadb a et aT~l ([2010), 
which is ba sed on the absorpt ion-corrected AGN X-ray 
luminosity (|Carilli et al. 112002 1 . Given the uncertainties 
in the IR and X-ray bolometric corrections, we consider 
the two estimates to be in reasonable agreement. Using 
our mid-IR based estimate of the hidden quasar luminos- 
ity and an assumed mass-energy conversion efficiency of 
0.1, we estimate a mass accretion rate of ~ 16M Q yr _1 
for the supermassive black hole that powers the Spider- 
web AGN. This is ~ 1% of the PAH-based instantaneous 
star formation rate integrated over the entire galaxy. 

4.2. HAE 229 and HAE 131 

HAE 229 (z = 2.149) and HAE 131 (z = 2.152) are two 
of the brightest Ha emission candi dates in the s urvey o f 
the PKS 1138-26 protocluster by iKurk et all (|2004aj) . 
and have confirmed Ha and [ N III emission in the ir NIR 
(rest- frame optical) spectra (|Kurk et al. I l2004b(l . The 
projected distances from the center of the Spiderweb 
Galaxy are 32'.'9 (280 kpc) and 49'.'9 (420 kpc), respec- 
tively. 

Note that the QSO CXO J114045.9-262916 (PKS 1138: 
[KPR2004] HAE 0215, z = 2.157) falls within the spec- 
tral extraction aperture 5"4 SW of HAE 229, (Fig. 1). 
The positional uncertainty for the QSO is 0'.'45, and we 
estimate an absolute positional uncertainty of < l'.'O for 
HAE 229 from our IRS spectral map. The two sources 
are spatially resolved at the wavelength of the 7.7 /.im 
PAH band. By comparing the PAH feature and con- 
tinuum extracted at the peak pixel of each source, we 
estimate that < 30% of the 7.7/zm PAH flux and < 60% 
of the 8 /im continuum in the HAE 229 spectrum can be 
contributed by the QSO. 

We measure a PAH luminosity for HAE 229 of 3.5 x 
lO 11 !/©. The inferred star-forming IR luminosity is 
Lsfir ~ 6.0 x 10 12 L Q and the PAH-derived star for- 
mation rate is ~ 880M Q yr _1 , similar to the Spider- 
web Galaxy. HAE 229 is one of two red galaxies in 
the PKS 1138-26 cluster that were detected in Hn by 
iDohertv etaTl (12010ft (galaxy 464). The Ha luminos- 
itv ((Kurk et al. ll2004bD indicates a star formation rate 
of = 20 Mq yr _1 , a factor of ~ 40 less than what we 
estimate from the 7.7 /im PAH feature (Table 3). This 
indicates that the bulk of star formation activity is heav- 
ily obscured (> 4.1 mag) at the wavelength of Ha. SED 
fitting yiel ds a stellar mass of 5 x 10 11 M Q , and an age 
of 2.4 Gyr ((Dohertv et al. 112010(1 . The time required to 
double the stellar mass is only 600 Myr, indicating that 
protocluster galaxy HAE 229 is undergoing a major event 
in its construction. 

We measure a PAH luminosity for HAE 131 of 1.9 x 
10 n i Q , half as bright as HAE 229 and the Spiderweb. 
The corresponding estimated IR luminosity and star for- 
mation rate are Lsfir ~ 3.2 x lO 12 ^ and ~ 470M© 
yr _1 . The Ha luminosity ((Kurk et al. ll2004bD indicates 
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a star formation rate of 34 M Q yr _1 , a factor of 14 less 
than what we estimate from the 7.7 /zm PAH feature 
(Table 3). As for HAE 229, most of the star formation 
in HAE 131 must be obscured (> 2.9 mag) at the wave- 
length of Ha. 

The ultraluminous PAH emission and large PAH/Ha 
flux ratios in both HAE 229 and HAE 131 clearly in- 
dicate that they are ULIRGs, undergoing a rapid burst 
of obscured star formation. Similar PAH luminosities 
and equivalent widths have been seen in ot her high- 
redshift, 24 /zm-selec ted ULIRGs (jSajina et al. I r 2007L 
iTanaka et all (|20 1 IT ) find that HAEs in the PKS 1138- 
26 cluster have systematically smaller Ha/K-band stel- 
lar continuum ratios than HAEs in the 4C 23.56 cluster 
(z = 2.483), attributing this to smaller specific star for- 
mation rates. However, we find Ha extinction factors 
in two PKS 1138-26 HAEs of 14 and 40, compared to a 
range of 2-14 in the 4C 23.56 cluster, based on Ha/24 /im 
fluxes. If this tendency towards larger extinction in PKS 
1138-26 HAEs holds up, the difference in Ha/K-band 
stellar continuum ratios may be ascribed to differences 
in extinction, rather than in specific star formation rates. 

The lack of high-ionization lines in the HAE 229 
spectrum indicates a relatively weak AGN contribution, 
though there is an indication of hot dust emission at at 6- 
7 /im rest. Detection of [Ne Vi] in the HAE 131 spectrum 
together with the detection of hot dust indicates an AGN 
contribution, though it must be relatively weak, consid- 
ering the large PAH equivalent widths. We do not detect 
H 2 0-0 S(3) emission in HAE 229 nor HAE 131, with a 2ct 
upper limit of < 1 x 10 -15 erg s _1 cm~ 2 (< 9 x 10 9 L & ), 
even though they have similar PAH luminosities and in- 
ferred star formation rates to the Spiderweb. This sup- 
ports the idea that the radio jet or AGN excites the H 2 
emission in the Spiderweb, rather than star formation. 



5. CONCLUSIONS 

We have found PAH emission indicating star forma- 
tion rates of 500 — 1100M Q yr _1 in three galaxies in 
the PKS 1138-26 radio galaxy protocluster, including 
the Spiderweb central radio galaxy, HAE 229, and HAE 
131. The corresponding stellar mass doubling time scales 
in the Spiderweb and HAE 229 are 1.1 Gyr and 600 
Myr, respectively. This indicates that these protoclus- 
ter galaxies are undergoing a period of rapid growth and 
that the PKS 1138-26 protocluster is a massive galaxy 
construction zone. Similarly high star formation rates 
have been observed in other high-redshift radio galaxies 
(e.g. iSevmour et al. 11200 8). in contrast to low-redshift 
[z < 0.2) radio galaxies where star formation rates are 
generally lower (< 5Mq yr _1 ), in spite of the presence 
of u p to 10 10 MfT) of molecular gas in the central few 
kpc llShi et al. 1(20071 lOele et al. Il20lot ISevmour et al. I 
1201 If ). Low-redshift radio jets appear to be more effec- 
tive at suppressing star-formation in their host galaxies. 
However, the dynamical state and gas physical properties 
of low-redshift radio galaxies are also much more evolved, 



as indicated by the presence of virialized hot gas in their 
halos and surrounding galaxy clusters. 

Star formation rates estimated from the 7.7 /zm PAH 
feature are a factor of 14-40 greater than estimates based 
on Ha, for HAE 229 and HAE 131. This indicates that 
star formation is heavily obscured (by > 2.9 — 4.1 mag) 
at Ha and shorter wavelengths. This is in contrast to the 
Spiderweb, where the PAH-based and Ha-based star for- 
mation rate estimates are comparable. This may suggest 
a different mode of more-extended, less obscured star 
formation in the Spiderweb. This in turn may be consis- 
tent with the observ ed extended UV continuum emission 
(jHatch et al. |[2008h if the UV extinction is more compli- 
cated than a simple obscuring screen. 

The extended mode of star-formation in the Spiderweb 
could possibly be a result of negative feedback quench- 
ing star-formation in the nucleus, coupled with accretion 
of many satellite galaxies enhancing star formation in 
the outer regions of the galaxy. High luminosity H 2 0-0 
S(3) emission from the Spiderweb indicates a large mass 
of warm (T > 300 K) molecular hydrogen, most likely 
heated by radio jet feedback from the AGN, which may 
have already quenched nuclear star formation and may 
eventually quench quasar activity. In other words, the 
high-star formation rates observed in high-redshift radio 
galaxies may be attributed to the accretion of gas and 
gas-rich satellite galaxies rather than nuclear star forma- 
tion. At low redshift, the gas has been virialized and 
heated to 10 7 K and the satellites have been incorpo- 
rated into the central radio galaxy, such that there is no 
longer a steady gas supply for star formation. The re- 
maining gas that is able to cool and make its way into the 
central kpc-scale disk is unable to form stars efficiently 
because it is either continuously or episodically heated 
by the radio jet. 

Hcrschcl FIR-sub-mm photometry will determine the 
total IR luminosities of the galaxies in the PKS 1138- 
26 protocluster, yielding more accurate star formation 
rates and dust mass estimates. Herschel PACS spec- 
troscopy may be sensitive enough to detect the H 2 S(0) 
line and thereby measure the total mass of warm molec- 
ular gas in the Spiderweb. The extreme luminosity of 
H 2 rotational line emission in the Spiderweb provides a 
strong incentive for launching a cooled, large aperture 
FIR space telescope to detect these lines and use them 
to measure the impact of kinetic energy dissipation on 
molecula r gas and star formatio n in even higher-redshift 
galaxies (jAppleton et al. 1 120091 ) . 

This work is based on observations made with the 
Spitzer Space Telescope, which is operated by the Jet 
Propulsion Laboratory, California Institute of Technol- 
ogy under NASA contract 1407. This work benefited 
greatly from the NASA/IPAC Extragalactic Database 
(NED), which is operated by the Jet Propulsion Labora- 
tory, California Institute of Technology, under contract 
with NASA. N. S. is a recipient of an Australian Research 
Council Future Fellowship. 
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TABLE 1 
PAH Features 



Source 



7.42 a 



7.60 



7.85 



8.33 



8.57 b 



11.23 



11.33 



Spiderweb 
HAE 229 
HAE 131 



2.30 (0.28) 

1.13 (0.57) 

< 0.56 



0.94 (0.08) 
0.95 (0.10) 



2.62 (0.07) 
1.88 (0.07) 
1.16 (0.07) 



0.28 (0.06) 
<0.14 



1.15 (0.04) 
0.62 (0.04) 
0.79 (0.06) 



1.08 (0.06) 
<0.14 
< 0.30 



0.29 (0.14) 
1.82 (0.16) 
2.24 (0.27) 



a Flux and la uncertainty or la upper limit (10 erg s cm — 10 W m ). 

PAH feature at 8.6 /Am is better fit using a central wavelength of 8.57 fim than the standard value of 8.61 ftm 









TABLE 2 
Emission Lines 








Source 


H 2 0-0 S(5)+[Ar n] 


[No vi] 


H 2 0-0 S(4) 


[Ar m] 


H 2 0-0 S(3) 


[Siv] 


Spiderweb 
HAE 229 
HAE 131 


0.80 (0.08) 
0.44 (0.07) 
0.17 (0.08) 


1.56 (0.04) 

< 0.30 
0.18 (0.04) 


0.61 (0.04) 

< 0.21 
0.20 (0.04) 


0.37 (0.03) 
< 0.07 


0.41(0.03) 

< 0.06 

< 0.11 


1.55 (0.04) 

< 0.17 

< 0.08 



Notes. Flux and la uncertainty or la upper limit (10 erg s cm — 10 W m ). 



TABLE 3 
Luminosities and SFRs 



Source 



L(PAH7.7) a 



L(Hcv) b 



SFR(PAH) d 



SFR(Ha) c 



Spiderweb 
HAE 229 
HAE 131 



4.4 (0.3) 


38 (3) 


7.4 (0.5) 


1100 (60) 


1200 (10 


3.5 (0.5) 


0.64 (0.17) 


6.0 (0.9) 


880 (130) 


20 (5) 


1.9 (0.2) 


1.1 (0.2) 


3.2 (0.3) 


470 (30) 


34(6) 



a .Luminosity (10 £©) of 7.7/im PAH blend, including contributions from the 7.42 |iin, 7.60 |im, and 7.85 /im features. 

b Ha luminosity (10 9 L Q ). Ha fluxes for HAE 229 and HAE 131 arc from IKurk et al. I ll2004bll . Spiderweb extended narrow Ha flux has AGN 

contribution subtracted {Ncsvadba ct al. 2006). 

c Estimated Star-forming total IR luminosity L S fir(10 12 L Q ) = 17 X L(PAH7.7) (Elbaz ct al. 2011). 

d PAH-based star formation rate SFR(PA H, M yr -1 )= L(PAH7 .7)/4 X 10 s (see text). 

° SFR(Ha, M yr _1 ) = 3.1 X 10~ 8 L(Ha) dKennicutt et al. |[l994T) . 

Figure 2 Caption 

Spitzer IRS LL1 spectra, (a) Spiderweb radio galaxy (z = 2.156). PAHFIT spectral fit (green) and fit components: 
continuum (solid black), silicate emission (dashed), PAH features (blue), and emission lines (purple), (b) Spiderweb 
PAH and emission line spectrum after subtracting continuum fit. Emission lines [Ar il]+H2 0-0 S(5), [Ne Vi], [Ar in], 
H 2 0-0 S(3), and [S iv]; and 7.7, 8.6, and the blue wing of the 11.3 ^m PAH are detected, (c) HAE 229 (z = 2.149) 
(d) HAE 131 (z = 2.152). PAH features are detected at 7.7, 8.6, and 11.3 ^m in both HAE 229 and HAE 131. The 
[Ar n] emission line is also detected. 
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Fig. 1. — Spitzer IRS maps (3.2' X 1.1') of the PKS 1138-26 protocluster. The mean slit orientation and long axis of each map is at 
PA= 125.1°. Top: IRS map in synthetic MIPS 24 (im band. The locations of spectroscopic cluster members (diamonds) or candidates 
(small circles) and X-ray sources (squares) are indicated. Spectral extraction apertures (large circles) are shown for the Spiderweb (center), 
HAE 229 (top), and HAE 131 (right). Note that the QSO CXO J114045.9-262916 (z = 2.157) falls inside and near the SW edge of the 
spectral extraction aperture for HAE 229, and may contribute flux to its spectrum (§4.2). Bottom: synthetic H2 0-0 S(3) narrow-band, 
continuum subtracted map. The same compass rose and scale apply to both maps. 
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Fig. 2. — See above for figure caption, which will not fit on the page. 
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Fig. 3. — Estimated mass of warm (300 < T < 1200 K) H2 in the Spiderweb Galaxy versus temperature (solid line). Temperatures 
T > 1200 K (dotted line) are excluded by the upper limit to th e H2 S(5) flux . I n comparison, warm molecular hydrogen masses in the 
same temperature range for z < 0.22 3C radio galaxy MOHEGs (Ogle et al. 2010) are also plotted (diamonds). 



